Out-of-equilibrium statistical ensemble inequivalence 
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We consider a paradigmatic model describing the one-dimensional motion of N rotators coupled 
through a mean-field interaction, and subject to the perturbation of an external magnetic field. 
The latter is shown to significantly alter the system behaviour, driving the emergence of ensemble 
inequivalence in the out-of-equilibrium phase, as signalled by a negative (microcanonical) magnetic 
susceptibility. The thermodynamic of the system is analytically discussed, building on a maximum 
entropy scheme justified from first principles. Simulations confirm the adequacy of the theoretical 
picture. Ensemble inequivalence is shown to rely on a peculiar phenomenon, different from the 
one observed in previous works. As a result, the existence of a convex intruder in the micro- 
canonical energy is found to be a necessary but not sufficient condition for inequivalence to be 
(macroscopically) observed. Negative temperature states are also found to occur. These intriguing 
phenomena reflect the non-Boltzmanian nature of the scrutinized problem and, as such, bear traits 
of universality that embraces equilibrium as well out-of-equilibrium regimes. 

PACS numbers: 



Classical statistical mechanics is most commonly deal- 
ing with large systems, in which the interaction range 
among components is much smaller than the system size. 
In such "short-range" systems, energy is normally addi- 
tive and statitistical ensembles are equivalent The 
situation may be radically different when the interac- 
tion potential decays so slowly that the force experienced 
by any system element is dominated by the interaction 
with far-away components. In these "long-range" inter- 
acting systems (LRI) energy is not additive. Well-known 
examples of non-additive LRI systems are for instance 
found in cosmology (self-gravitating systems) and plasma 
physics applications, where Coulomb interaction are at 
play 0, [J- The lack of additivity, together with the pos- 
sible non-convexity of the space of accessible macroscopic 
thermodynamic parameters (break of ergodicity) , may be 
at the origin of a large gallery of peculiar thermodynamic 
behaviours: the specific heat can be negative in the mi- 
crocanonical ensemble 0, and temperature jumps may 
appear at microcanonical first-order phase transitions. 
When this occurs, experiments realized on isolated sys- 
tems give different result from similar experiments per- 
formed on systems in contact with a thermal bath. As 
a consequence, the canonical and microcanonical statis- 
tical ensembles of long-range interacting systems may be 
non-equivalent . 

Ensemble inequivalence has been theoretically ob- 
served in several models of LRI systems, at thermody- 
namical equilibrium 0. In all cases, for isolated systems 
such an observation has been related to the a convex 
shape of the microcanonical entropy Indeed, if the 
entropy s(e) is twice differentiable and non-concave for 
some values of the specific energy e, the microcanonical 
specific heat is negative. Thus, since the canonical (ther- 
modynamic) specific heat is a positive-defined quantity, 
ensembles are non-equivalent. In open systems, s may 



also depend on other extensive thermodynamic variables, 
apart from energy. In this case, (equilibrium) ensemble 
inequivalence has been detected also if s is non-convex 
along e, provided it is convex along one of the other vari- 
ables [f|. 

Systems with LRI also display a quite peculiar dynam- 
ics, which is characterized by a slow relaxation towards 
thermodynamical equilibrium, and, even more remark- 
ably, by the convergence to out-of-equilibrium quasi- 
stationary states (QSS) It has been shown that 
QSS's can be related to the stable steady states of the 
Vlasov equation, which describes the system in the limit 
TV — > oo [3|. The idea was inspired by the seminal work 
of Lynden-Bcll H , developed in the context of stellar dy- 
namics, and later applied to vortex dynamics (see, e.g., 
) ■ As it will be discussed more in detail in the follow- 
ing, Lynden-Bell's approach is based on the definition of 
a locally-averaged ("coarse-grained") distribution func- 
tion, which then translates into an entropy functional, as 
follows from standard statistical mechanics prescriptions. 
By maximizing such an entropy, while imposing the con- 
straints of the dynamics, returns a closed analytical ex- 
pression for the single particle distribution of the system 
in its QSS regime. The microcanonical Lynden-Bell en- 
tropy enables one to introduce an out-of-equilibrium free- 
energy, defined via the usual Legendre-Frenchel trans- 
form [H, [Til] - This makes it possible to estimate the 
relevant collective variables, as calculated both in the 
microcanonical and canonical out-of-equilibrium statisti- 
cal ensembles, to eventually challenge the equivalence of 
the two formulations. Following this procedure, it was 
recently demonstrated [Tl| that (out-of-equilibrium) en- 
semble equivalence may still hold when negative kinetic 
specific heat is measured in both statistical ensembles. 
While QSS's represent out-of-equilibrium states of the 
Af-body dynamics, they could be equally interpreted as 
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equilibrium configuration of the corresponding continu- 
ous description: in this respect, the conclusions of our 
analysis will apply to both equilibrium and non equi- 
librium dynamics, provided the latter bears distinctive 
non-Boltzmannian traits. 

In this letter, we extend the concept of statistical in- 
equivalence to out-of equilibrium LRI systems. As it will 
be shown, the inequivalence does not materialize as an 
immediate byproduct of a "convex intruder" in the micro- 
canonical entropy. Indeed, the entropy is always a con- 
vex function of one of the thermodynamic variables, over 
the whole domain of its definition, including the region 
where experiments realized in the two ensembles would 
provide the same results (ensemble equivalence). This 
observation, that we here substantiate via numerical and 
analytical means, suggests that the presence of a "convex 
intruder" in the entropy is a necessary but not sufficient 
condition for the inequivalence of statistical ensembles to 
be observed in a macroscopically tangible form, e.g. as 
a difference in sign of the system's specific heat and/or 
susceptibility. This phenomenon should be therefore cat- 
egorized as a new type of inequivalence, distinct from 
the one so far reported in the literature. Further, the 
inequivalence is here shown to emerge as an effect of an 
externally imposed forcing, that significantly impacts the 
associated, unperturbed dynamics for which equivalence 
is shown to hold Finally, the results here derived 

can be also cast in the framework of an equilibrium non- 



Boltzmanian picture, as in the spirit of [11 



As a paradigmatic example for our investigation, we 
consider a model describing the one-dimensional motion 
of N rotators coupled through a mean-field cosine in- 
teraction and subject to the perturbation of an external 
magnetic-like field. The model is mathematically defined 
by the following Hamiltonian: 
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(1) 

Here 9j represents the orientation of the j-th rotor 
and pj is its conjugate momentum. The parameter h 
is a scalar quantity, which measures the strength of the 
magnetic field. For h = 0, the previous Hamiltonian re- 
duces to the celebrated Hamiltonian Mean Field (HMF) 
model [l2l |. which has been widely studied in the past 
as a prototype of LRI system. To monitor the evo- 
lution of the system, it is customary to introduce the 
magnetization, a macroscopic order parameter defined 
as M = |M| = |^nij|/./V, where m.; = (cos0j,sin#j) 
stands for the microscopic magnetization vector. 

As previously reported [13], for h = 0, after an initial 
transient, the system gets trapped into QSS's. Similarly, 
QSS's are also encountered in the generalized setting with 
h ^ 0. 

In the N — >• oo limit, i.e. when the system is indefi- 
nitely stuck in the QSS's, the discrete dynamics gener- 
ated by the Hamiltonian ((TJ) can be described in terms 



of a continuum Vlasov equation. The QSS's are then in- 
terpreted as the stable (attractive) steady states of the 
underlying Vlasov equation. Analytical progress is possi- 
ble by invoking the aforementioned Lynden-Bell violent 
relaxation theory Q, that we shall tackle with reference 
to a simplified choice of the initial condition. Assume 
that the initial single particle distribution takes only two 
distinct values, namely fo and zero. Angles and veloci- 
ties populate a bound domain in phase space, therein dis- 
playing a uniform probability distribution. This working 
ansatz corresponds to dealing with the so-called "water- 
bag" distribution. Vlasov time evolution can modify the 
shape of the boundary of the "water-bag", while con- 
serving the area inside it. Thus, the distribution remains 
two-level (0, /o) as time progresses. By performing a lo- 
cal average of / inside a given mesoscopic box, one gets a 
coarse-grained distribution /, which eventually converges 
to an asymptotic equilibrium profile. The latter solution 
is explicitly calculated by maximising an entropy func- 
tional, associated to /, through a direct combinatorial 
calculation that moves from a rigorous statistical me- 
chanics setting Q. In the two- level scenario, the mixing 
entropy s can be cast in the form: 
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If one operates in the setting of Hamiltonian (fTJ, 
both the specific energy e(/) = / j (p 2 /2)f(9,p,t)d9dp - 
{Ml + M 2 - l)/2 - hM x and the momentum P(f) 
j J pf(9,p, t)d9dp are conserved by the dynamics. In ad- 
dition, the normalization of the distribution / has to be 
accounted for, which in turn corresponds to dealing with 
constant mass. Requiring the entropy to be stationary, 
while imposing the conservation of the above quantities 
(e(/) = e and P{f) = cr), defines a variational problem 
that admits as a solution the following fermionic distri- 
bution: 



/qss(0,p) 



1 + e /3/o(p 2 /2 - M[/ QSS ] • m - hcos8) + a ' 



(3) 



where the label QSS recalls that the recovered Vlasov 
equilibrium distribution is meant to describe the out-of- 
equilibrium (QSS) iV-body regime. Here a and /3 play 
the role of Lagrange multipliers associated, respectively, 
to mass and energy conservation; M = (M x ,M y ), where 
M x [f] = J fcos(9)d9dp, M y [f] = J fsm{6)d6dp stand 
for the two components of magnetization in the N — > oo 
limit; finally, m = (cos 9, sin 9). Once the system energy 
e and fo are fixed, one can determine the correspond- 
ing values of M x , M y , a and /3 in the expression ((3|), 
by self-consistently imposing the constraints condition. 
The obtained system of equations is then solved using a 
Newton-Raphson method. 

We are in particular interested in monitoring the re- 
sponse of the system to the external solicitation here 
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FIG. 1: Magnetization M in the QSS, calculated using the 
Lynden-Bell approach, plotted as a function of h, for different 
choices of e (continuous and segmented lines). From top to 
bottom: e = 0.51,0.53,0.62,0.64. In the inset, the results 
of direct iV-body simulations (symbols) are compared to the 
Lynden-Bell prediction, for e = 0.53 and e = 0.62. In the 
simulations, N = 10000. Each data point refers to an average 
over 20 independent realizations. The error bar is calculated 
as the variance of the 20 independent estimates. 



encoded in the parameter h. Fixing the energy, while 
changing h returns a sequence of distinct equilibrium 
values of the magnetization, as calculated via the above 
maximum entropy procedure. It is worth noting that the 
stationary entropy depends paramctrically on both the 
energy and (quasi-stationary) magnetization, the latter 
being independent thermodynamic quantities. The mag- 
netization values are displayed in fig|T]as a function of h, 
and for different energy choices. It is worth emphasising 
that we here operate at constant / . Surprisingly, below 
a threshold in energy the system shows negative suscepti- 
bility x = 9M /dh, the magnetization decreasing with the 
strength of the externally imposed field. Above the crit- 
ical energy value, the magnetization grows instead with 
h, which implies dealing with a positive susceptibility. 

The presence of a region with negative susceptibil- 
ity in the microcanonical system (JlJ, points to out- 
of-equilibrium ensemble inequivalence. Imagine in fact 
to introduce the system free energy, by performing the 
Legendre-Frenchel transform of entropy ((2]) [ll| . In the 
so-defined canonical ensemble, the sign of x is related to 
the concavity of the free energy. The latter being (by 
definition) always concave, the canonical susceptibility is 
positive defined. In other words, x can be negative in the 
microcanonical ensemble, as it does in a specific param- 
eters region, while it is always positive definite quantity 
in the canonical framework. The theoretical evidence for 
X < represents hence an indirect (macroscopically ac- 
cessible) signature of ensemble inequivalence. 

To confirm the adequacy of the proposed theoretical 
picture we performed a campaign of direct N-body simu- 
lations, sampling the QSS, after the initial transient has 
died out. The magnetization is computed as a function 
of h, averaging over many independent realizations. The 
recorded data are plotted in the inset of figffl showing a 
quite satisfactory agreement with the analytical calcula- 
tion. 



In the remaining part of the paper, we shall elaborate 
on a thermodynamical description of the phenomenon 
here detected, by exploring the topological properties of 
the entropy functional ^ versus energy and magnetiza- 
tion. 

For a magnetic system, the the first principle of ther 
modynamics can be cast in the form 
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Tds = de + hdM 



(4) 



where T is the thermodynamic temperature and e is 
the total energy of the system, including the field contri- 
bution . Recalling that s is a function of the independent 
variables (e, M) one readily obtains the following expres- 
sion for h: 



V dM I e 
\deJ M 



(5) 



where use has been made of the relation T = 
(ds/de)^ . From the expression for h, one can imme- 
diately cast the susceptibility in the micro-canonical sta- 
tistical ensemble as: 
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where for the sake of clarity we dropped the label that 
identifies the variable kept constant upon differentiation. 
The sign of x, and thus the associated (macroscopically 
accessible) out-of-equilibrium inequivalence features, are 
ultimately controlled by the (first, second and crossed) 
partial derivatives of the entropy s with respect to M 
and e. Using the Lynden-Bell approach, the three di- 
mensional profile of s can be analytically reconstructed 
as a function of the energy e and (sclf-consistcntly de- 
termined) magnetization M. In fig. [5] we report a 2D 
projection of the entropy iso-h lines in the (M, e) plan, 
for fo = 0.12: each line refer to a specific value of h, 
which is frozen to the selected value, while visiting the 
entropy surface. The lines cross in a well defined point of 
the plan that we identify as (M c , e c ) ~ (0.46, 0.575) fl5| . 
Two different domains can be singled out, which, as will 
become transparent later on, appear to be disconnected: 
a larger region, termed I in the following, delimited by 
the conditions e > e c and M < M c . A smaller region II, 
in which e < e c and M > M c [lj| . 

In figOk, the entropy s is plotted versus M while 
keeping e constant: a clear qualitative change is de- 
tected when going from region I (dashed lines) to re- 
gion II (solid lines). In the latter domain, the entropy 
decreases with the magnetization M, while in the for- 
mer the entropy grows as M gets larger. In region I, 
hence (ds) / '(dM) > 0: from relation ([S]), and having 
chosen positive values of h, we expect (ds)/(de) > 0, 
which implies working with positive temperatures. This 
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FIG. 2: Projection of s(e,M) in the (M, e) plan: each line 
is an iso-h, with h — 0.1,0.2,0.3,0.4, from bottom to top 
in region I (vice versa in region II). The lines intersect in 
a unique point, (M c ,e c ) ~ (0.46,0.575). Here, f = 0.12. 
Regions I and II as defined in the text are visually delimited 
by two solid, orthogonally crossing, lines. 




What can we say about the susceptibility \ l The lat- 
ter can be calculated using relation ©. The result of the 
calculation are reported in figure^] where the outcome of 
the analysis in a sub-domain of region I and II is plotted 
in separate panels. As confirmed by direct inspection, in 
region I the susceptibility is positive, while it is negative 
in region II. In the latter domain, bound in energy from 
below, the system thus displays ensemble inequivalence, 
as revealed by the change in sign of a macroscopically 
accessible thermodynamic quantity. Intriguingly enough, 
inequivalence is also accompanied by the additional un- 
conventional condition T < 0. Coming back to figO it is 
worth emphasising that the entropy is a concave function 
of the energy e (at constant M), while it is always non 





FIG. 4: The susceptibility \ is plotted as a function of both 
M and e, in region I (fig. a)) and region II (fig. b)). 



FIG. 3: Fig. a): Entropy s is plotted versus (QSS) mag- 
netization M, for different choices of the energy: e = 
0.51,0.53,0.62,0.64, from left to right. Fig. b): s is plot- 
ted versus e, for different choices of the (QSS) magnetization 
M: M = 0.49, 0.51, 0.45, 0.43, from left to right. 

prediction is confirmed by visual inspection of figOjD: in- 
deed, the entropy s increases with the energy e, when 
sampling the portion of the plan associated to region I. 
More interesting it is instead the scenario that emerges 
from a close look of region II. Here, (ds)/{dM) < and 
T _1 = (ds)/(de) < 0. Summing up the information 
condensed in figure [3J the temperatures in region I are 
positive, while those associated to region II are nega- 
tive. When the critical point (M c , e c ) is approached the 
entropy diverges, an observation that suggests targeting 
region I and II as disconnected domains. 



concave as a function of the magnetization M (at fixed 
e): as we have demonstrated, such a feature is respon- 
sible for ensemble inequivalence in region II, but is not 
sufficient to induce observable inequivalence in region I. 

In this letter we have studied a paradigmatic model 
of LRI system, perturbed by an external magnetic field. 
By using a maximum entropy technique grounded on first 
principles, we have given both analytical and numerical 
evidence of (out-of-equilibrium) negative susceptibility. 
Such a result points to the existence of out-of-equilibrium 
ensemble inequivalence. According to our analysis, and 
at difference from previously reported results, the exis- 
tence of a convex intruder in the entropy does not trans- 
late in a sufficient criterion for the negative susceptibility 
to manifest. 
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